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Abstract 



We illustrate the sensitivities of LEP experiments to leptonic signals 
associated to models where supersymmetry (SUSY) is realized with spon- 
taneous breaking of R-parity. We focus on missing transverse momentum 
plus acoplanar muon events {'pT + A*^/^ ) arising from lightest neutralino 
single production well as pair production xx, followed by x decays, 

where x denotes the lightest neutralino. We show that the integrated lu- 
minosity achieved at LEP already starts probing the basic parameters of 
the theory. We discuss the significance of these constraints for the sim- 
plest spontaneous R-parity breaking models and their relevance for future 
searches of SUSY particles. 



*E-mail fromao@alfa.ist.utl.pt 
t E-mail fernando® axp . ift . unesp . br 
■l-E-mail miguel@flamcnco.ific.uv.es 
^E-mail magro@flamenco.ific.uv.es 
^E-mail valle@flamenco.ific.uv.es 



1 Introduction 



So far most searches for supersymmetric particles have been made in the framework of 
the Minimal Supersymmetric Standard Model (MSSM) which assumes the conservation 
of a discrete symmetry called R-parity [Q]. Under this symmetry all the standard model 
particles are R-even while their superpartners are R-odd. R-parity is related to the spin 
(S), total lepton number (L), and baryon number (B) according to Rp = (— 
In the limit of exact R-parity the supersymmetric (SUSY) particles must be produced 
only in pairs, the lightest of them being stable. 

Unfortunately there is no clear dynamical clue as to how supersymmetry is re- 
alized. In fact, neither gauge invariance nor supersymmetry require the conservation 
of R-parity. The violation of R-parity could emerge as the residual effect of a more 
fundamental unified theory . While R-parity violation may be explicit , we find it 
rather attractive to consider the possibility that it arise in a spontaneous way, like the 
breaking of the electroweak symmetry [^]. At the present state-of-the-art theory can 
not decide. It is therefore of great interest to pursue the phenomenological implica- 
tions of alternative scenarios. This is specially so in view of the fact that the associated 
phenomena can be accessible to experimental verification. 

Recently there was a lot of attention devoted to the possibility that R-parity 
can be an exact symmetry of the Lagrangian, broken spontaneously through nonzero 
vacuum expectation values (VEVS) for scalar neutrinos H, 0, ||, |^, |10|, |Tl], [1^. There 



are two main types of scenario |jio[. If lepton number is part of the gauge symmetry 
there is an additional gauge boson which acquires mass via the Higgs mechanism. In 
this case there is no physical Goldstone boson and the scale of R-parity violation also 
characterises the new gauge interaction, around the TeV scale Consequently, its 
effects can be large 0. In this model typically the lightest SUSY particle (LSP) is a 
neutralino which decays mostly to visible states, breaking R-parity. The main decay 
modes are three-body, 

X - //^ (1) 
where f denotes a charged fermion. Its invisible decay modes are in the channel 

X^3u (2) 

Alternatively, if spontaneous R-parity violation occurs in the absence of any additional 
gauge symmetry, it leads to the existence of a physical massless Nambu-Goldstone 
boson, called Majoron (J) |T^. Thus in this case the lightest SUSY particle is the Ma- 
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joron which is massless and therefore stable [jj]. As a result often the lightest neutralino 
X niay decay invisibly, conserving R-parity, as 

X^i^ + J. (3) 

With the minimal particle content of the MSSM, including only the usual isodoublet 
sneutrinos, these models lead to a new decay mode for the Z boson, Z p + J, where 
p is a light scalar with mass -C M^. This decay increases the invisible width of the Z 
by the equivalent of one half extra neutrino species, and is therefore ruled out by the 
LEP measurements. This difficulty, as well as the fine-tuning problem characteristic of 
this simplest scenario, is naturally avoided by adding SU{2) ® f/(l) singlet neutrinos, 
in such a way so that the R-parity breaking is driven by the corresponding sneutrino 
VEVS. In this case the Majoron is mainly singlet |[. 

In this paper we analyse some aspects of the phenomenology of spontaneously 
broken R-parity models at the Z peak. We consider the production of the lightest su- 
persymmetric fermion, including both single as well as pair production mechanisms. In 
order to identify the corresponding signals we take into account all possible neutralino 
decay channels, thus extending previous discussions. We focus on the study of events 
with missing transverse momentum plus acoplanar muons {"^t + f^^f^~), arising both 
from neutralino single production at the Z peak 

Z^X^ , (4) 
as well as pair production mechanisms 

Z ^XX, (5) 

followed by the R-parity violating x decays in eq. (p or eq. (Q) and eq. (^. Using 
the integrated luminosities attained by the four LEP experiments we conclude that 
the basic parameters of the model are starting to be probed in a meaningful way. 
We discuss the theoretical significance of these constraints from the point of view of 
other R-parity-violating processes as well as their possible implications for future SUSY 
particle searches. 



"The majoron may have a small mass due to explicit breaking effects at the Planck scale. In this 
case it may decay to neutrinos and photons. However the time scales are only of cosmological interest 
and do not change the signal expected at the laboratory [O. 
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2 Basic Structure 



Most of our subsequent analysis will be very general and applies to a wide class of 
SU{2) (g) f/(l) models with spontaneously broken R-parity, such as those of ref. |[, 
as well as models where the majoron is absent due to an enlarged gauge structure 
P, 0. Many of the phenomenological features relevant for the LEP studies discussed 
here already emerge in an effective model where the violation of R-parity is introduced 
explicitly through a bilinear superpotential term of the type iHu |1T6|. 



For concreteness, we start by adopting the conceptually simplest model for the 
spontaneous violation of R-parity proposed in ref. |^ and by recalling its basic ingre- 
dients. This will serve to set up our notation in what follows. We consider this model 
as the most useful way to parametrise the resulting physics, due to the strict correla- 
tion exhibited in this model between the magnitude of R-parity violating phenomena 
and the u^- mass. For example, all single SUSY particle production rates as well as the 
lightest neutralino decay rate are directly correlated to the mass of the tau neutrino. 
The superpotential is given by 

huQHuu" + hdHdQd" + hdHdc" + {hoH^Hd - e')$ + h^Huy" + h^Su' + h.c. (6) 

This superpotential conserves total lepton number and R-parity. The superfields ($, 
v'^i, Si) are singlets under SU{2) U{1) and carry a conserved lepton number assigned 
as (0, —1, 1) respectively. The couplings h^, h^, h^, h are arbitrary matrices in gen- 
eration space. The additional singlets v'^, S and $ [|1^] may drive the spontaneous 



violation of R-parity. This leads to the existence of a Majoron given by the imaginary 
part of 

'VuHu - VdHd) + —Ur - -TT^ T + -TjSr (7) 



Vv'^^ ' V V V 

where the isosinglet VEVS 



with V = yvj^ + Vg characterise R-parity or lepton number breaking and the isodoublet 
VEVS 

Vu = {Hu) , Vd = (Hd) (9) 

are responsible for the breaking of the electroweak symmetry and the generation of 
fermion masses. The combination v'^ = v^^ + v1 is fixed by the VF, Z masses. Finally, 
there is a small seed of R-parity breaking in the doublet sector, i.e., 

VL = {i>Lr) (10) 
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whose magnitude is now related to the Yukawa couphng hi^. Since this vanishes as 
^ 0, we can naturally obey the limits from stellar energy loss ||19|. 



For our subsequent discussion we need the chargino and neutralino mass matrices. 
The form of the chargino mass matrix is common to a wide class of SU (2) ® f/(l) SUSY 
models with spontaneously broken R-parity and is given by 









-iW+ 




heijVd 


— KijVRj 


V2g2VLi 


Hd 


heijVLi 


fi 


V^g2Vd 


-iW- 





V^g2Vu 


M2 



'111 



Two matrices U and V are needed to diagonalise the 5x5 (non-symmetric) chargino 
mass matrix 



xt = v.ji^t 
Xi = Uiji)j 



(12) 
(13) 



where the indices i and j run from 1 to 5 and tp'^ = {ef ,62 , H+, —iW'^) and 
V'j" = {eI,e2,e3,H^,-iW-). 

Under reasonable approximations, we can truncate the neutralino mass matrix 
so as to obtain an effective 7x7 matrix of the following form H 





Vi 


Hu 


Hd 




-iB 


J^i 





KijVRj 





92VLi 


-QlVLi 


Hu 


huijVRj 





-/i 


-g2Vu 


giVu 


Hd 





-/^ 





g2Vd 


-giVd 




92VLi 


-92Vu 


92Vd 


M2 





-iB 


-giVLi 


giVu 


-giVd 





Ml 



(14) 



This matrix is diagonalised by a 7 x 7 unitary matrix N, 



(15) 



where ifj^j = {ui, Hu, Hd, —iWs, —iB), with z/j denoting the three weak-eigenstate neu- 
trinos. 

In the above two equations Mi 2 denote the supersymmetry breaking gaugino 
mass parameters and gi^2 are the SU{2) ® U{1) gauge couplings divided by ^/2. We 
assume the canonical relation M1/M2 = ^tan'^Ow- Note that the effective Higgsino 
mixing parameter /i may be given in some models as fi = K {^), where ($) is the VEV 
of an appropriate singlet scalar. 
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Typical values for the SUSY parameters /i, M2, and the parameters /ij^j^s lie in 
the range given by 



(16) 



-250 < < 250 30 < < 1000 

10-10 < Kis, K23 < 10-1 10-5 < /i,33 < 10-1 

while the expectation values lie in the range: 

vl = VL3 = 100 MeV VLi = vl2 = 

50 GeV <vr = vr3 < 1000 GeV vri = vr2 = (17) 

50 GeV <vs = vs3 = vr < 1000 GeV 1 < tan /3 = ^ < 30 

The diagonalization of eq. (|14D gives rise to the mixing of the neutralinos with the 
neutrinos, leading to R-parity violating gauge couplings. In what follows we will give 
explicit expressions for the couplings of the SUSY fermions in terms of these diagonal- 
izing matrices. 



3 Charged and Neutral Current Couplings 

Using the above diagonalizing matrices U, V and N one can write the electroweak 
currents of the mass-eigenstate fermions. For example, the charged current Lagrangian 
describing the weak interaction between charged lepton/chargino and neutrino/neutra- 
linos may be written as 

-^W^XTl^iKLikPL + KRikPR)xl + H.C. (18) 

where Pl,r are the two chiral projectors and the 5x7 coupling matrices K^ r may be 
written as 

3 

KLik = m{-^Ui^NkQ - - E U^mNkm) (19) 

m=l 

KR,k = ek{-V2Vi,Nke + VMi) (20) 
The matrix -ftTiifc is the analogous of the matrix K introduced in ref |2^ and there 



is a corresponding matrix characterizing right-handed charged currents KRik. These 
"off-diagonal" blocks corresponding to i = 1..3 and k = 4. .7 as well as i = 4, 5 and 
k = 1..3 are R-parity-breaking couplings. 



The corresponding neutral current Lagrangian may be written as 

-z,{x7Yio'L,,PL + o'^,^Pr)x-u + \xh'{0'i,PL + o'^,,^Pr)xI} (21) 



COS Uy^r Z 



where the 7x7 couphng matrices p, and 0'[ ^ are given by 

O'lik = ViVk (^UiiUki + Ui5Uk5 + \Y1 UimUkm " ^ifc sin^ 6'vi/ j 



bv 




5ik sin^ Ow j 


(22) 


- sin^ 9w 


(23) 




(24) 



The off-diagonal part of these couphng matrices break R-parity, i.e. when i = 4.. 7 and 
= 1..3 or vice- versa. 

In writing these couphngs we have assumed CP conservation. Under this assump- 
tion the diagonahzing matrices can be chosen to be reaL The rji and factors are sign 
factors, related with the relative CP parities of these fermions, that follow from the 
diagonalization of their mass matrices. 

Like all supersymmetric extensions of the standard model, the spontaneously 
broken R-parity models are constrained by data that follow from the negative searches 
for supersymmetric particles at LEP, in particular the most recent limits on chargino 



masses from the recent run at 130 GeV mass region as well as pp collider data 
gluino production ^2 . 



There are additional restrictions, which are more characteristic of broken R-parity 
models. They follow from laboratory experiments related to neutrino physics and 
weak interactions, cosmology and astrophysics [p^ p!3|. These restrictions play a very 



important role, as they exclude many parameter choices that are otherwise allowed by 
the collider constraints, while the converse is not true. The most relevant constraints 
come from neutrinoless double beta decay and neutrino oscillation searches, direct 
searches for anomalous peaks at vr and K meson decays, the limit on the tau neutrino 
mass cosmological limits on the Ur lifetime and mass, as well as limits on muon 



and tau lifetimes, on lepton flavour violating decays and universality violation. 

One can perform a sampling of the points in our model which are allowed by all 
of the above constraints in order to evaluate systematically the attainable value of the 



couplings |2|]. The allowed values for the diagonal (R-parity conserving) couplings for 
the lightest neutralino x ^i-nd the lightest chargino ^ire of the same order as those 
in the MSSM. 

For the neutral current couplings of the lightest neutralino x we can only get an 
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upper limit After imposing the experimental constraints explained before, we get 



10^441^0.1 (25) 

In what concerns the R-parity breaking couplings, the biggest ones correspond to the 
standard lepton belonging to the third family, i.e. 0^43, and this will be responsible 
for the Z xv^ decay. This coupling amplitude can reach a few per cent or so for 
neutralino mass values accessible at LEP and, as we will show, may lead decay 
rates observable at LEP. 

In what follows we will focus on some of the zen-event signals that could be 
associated to neutralino single as well as pair production and subsequent decays at 
LEP. 



4 Neutralino Production at the Z Peak 



At the LEP I collider the neutralinos may be produced as 



e e 



(26) 



The differential cross section for these processes including only the Z exchange at 
-^i = Mz is given by 



do 



,e e 



XiXj) 



4J2 



i2-S.,)\Qis)\'\'/' I 



1 1 
s s 



sin 6^, cos 9^1 



Giij{s) + G2ij{s) cos 9 + Gsij{s) cos^ 9 



where A is the usual Kallen function and 

Ei E. 



Giij{s) 

G2ij{s) 

Gzijis) 



9v + 9l 



29V9A {02, - 02j) A^/^ (1 

l{9'v + 9\) {o2, + o2,)x{i 



s s 



and 



Q{s] 



s + mj — rn?- 
, Ej 



' 1 1 
s s 



s + — mf 



(27) 

(28) 
(29) 
(30) 



2v^ 



2v^ 



s - M| + iMzTz 



(31) 
(32) 
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Here gv and gA are the usual vector and axial couplings for the Standard Model Ze'^e 
vertex 

9v = + sin^ 9^ ; g^ = -- (33) 

and the relevant Xi Xj coupling amplitudes are determined from eq. (^). Notice that, 
due to the Majorana nature of the neutralinos, their neutral current couplings obey 
the relation O'^^j = —tiejO'^^j and therefore the functions G2ij{s) vanish identically. In 
the above equations the indices i and j run from 1 to 7. At = Mz they will be 
restricted only to those corresponding to neutralinos lighter than the Z boson. 

As seen above, in models with spontaneously broken R-parity, the mixing of 
the standard leptons with the supersymmetric charginos and neutralinos leads to the 
existence of R- violating couplings in the Lagrangian when written in terms of the mass 
eigenstates. As a result, SUSY particles can be singly-produced. This means that 
Xi and Xj sq. (^61) can be both supersymmetric particles (standard SUSY pair 
production) as well as one standard and one supersymmetric (R-parity breaking single 
production). This is in sharp contrast with explicitly broken R-parity models such as 
considered in ref. , where only the pair production of the lightest neutralino x can 
take place at the Z peak. 

In this paper we are therefore concerned with the signals arising from Z — 5> ^z/ 
decays (corresponding to i = j = 4) and Z — > XX decays (corresponding to i = 4, j = 
3 or vice- versa), where x denotes the lightest neutralino. The heavier ones will be 
assumed to be too heavy to be produced at the Z peak. 

5 Neutralino Decays 

Once produced in e"'"e~ collisions, the neutrahno subsequently decays, typically inside 

the detector. In order to identify the expected signatures at LEP it is necessary to 

specify its possible decay modes. In the MSSM all supersymmetric particles have 

cascade decays finishing in the LSP which is normally a neutralino. However, if R- 

parity is broken there are new decay channels and the supersymmetric particles can 

decay directly to the standard states breaking R-parity. Also the lightest SUSY particle 

may not be a neutralino, in fact this is the case in spontaneously broken models |^ 

Alternatively, SUSY particles may decay through R-parity conserving cascade decays 

** In this case the LSP may be produced in R-parity-violating decays of normal particles like the 
muon, the tau lepton jl^ or the Z boson pO| . 
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that will produce the lightest neutralino x- This may decay invisibly conserving R- 
parity as in eq. or via the R-parity breaking three fermion modes in eq. (|l]). 

To the extent that the invisible dominates, as could happen, one recovers the lightest 
neutralino missing momentum signal expected in the MSSM. 

For simplicity we are going to study the decays of the lightest neutralino, which 
one expects would be the earliest-produced supersymmetric particle. Heavier states 
would have cascade decays that we are not going to consider here. 

The lightest neutralino x has to decay always to standard states breaking R- 
parity. If its mass is lower than the mass of the gauge bosons it decays to the three 
body final states. We will distinguish three cases. First all the particles in the final 
state are neutral. Then 

X° ^ Vj Vk Vn, with width V\ = T^^{M^o, O'i O'^) (34) 

In the second case we consider charged leptons in the final state, that is 

/ ^ I- 1+ with width = T^''{M^o, 01 01 O'l, O'ji, Kl, Kr) (35) 

In both these processes there is interference among the various diagrams contributing. 
This is clear for the process in Eq. ^ where there is interference between the neutral 
and charged current diagrams, but it also true for the process in Eq. ^ due to the 
Majorana nature of the neutrinos. The explicit expressions for the widths r^** and F'^'" 
are given in the Appendix. Finally there is a third case when the final state contains 
quarks, that is 

X° ^i^.qq with width F" = V^'\M^o, O'i O^ O^ O'^, 0, 0) (36) 

that proceeds only via neutral current. In the Appendix it is explained how this width 
can be obtained particular case of F^ . 

As mentioned in the introduction, the existence of the Majoron implies that in 
SU{2)^U{1) spontaneously broken R-parity, the neutralino can always decay invisibly 
eq. d^) with a decay width 

CL^J = -e.e,Cffi, = pe,{N,kN,, + N,,N,,)Kk3^ ^^^^ 

Although our discussion will be more general, we neglect, for definiteness, supersym- 
metric fermion decays mediated by slepton exchange. In this approximation, neutrali- 
nos of mass accessible at LEP have only three-body decay modes mediated by charged 
and neutral currents, except for the two-body majoron decay eq. (|]), characteristic of 
the simplest spontaneous R-parity breaking models. 



9 



6 Signals at LEP 



In order to study the experimental signals associated to the first kinematically ac- 
cessible neutral supersymmetric fermions, we have developed an event generator that 
simulates the processes expected for the LEP collider at ^/s = Mz- It allows us to 
estimate the detection efficiencies when suitable selection criteria are imposed in order 
to avoid the expected standard model backgrounds to the processes of interest. We 
describe below the main steps we follow in order to generate neutralino production and 
decays. As far as the production is concerned, our generator simulates the following 
processes at the Z peak: 

• a) e~^e~ — *■ 

• b) e+e" XX 

Process a) clearly violates R-parity, so it is a new mode of neutralino production 
forbidden in the MSSM, as well as in models of explicitly broken R-parity other 
than those where this violation is due to a bilinear superpotential term iHu ll^. On 
the other hand, process b) is allowed both in the MSSM as well as in models such as 
the one used in ref. 

The second step of the generation is the decay of the lightest neutralino, which 
is the most characteristic feature of the R-parity breaking models. As explained in the 
introduction, if the lightest neutralino x is lighter than the Z boson it will have three 
body decays via charged or neutral currents, as well as the two body invisible decay 
into neutrino + majoron. These decays produce new supersymmetric signals and the 
generator allows their detailed study. In contrast, the two-body neutralino decay into 
neutrino + majoron has as signal missing transverse momentum, because both final 
particles escape detection and therefore we do not need to generate this process. Thus 
it suffices for us to generate the three body neutralino decays: 

• X — s> Ut-Z* Z/t-Z/I/, frQi^li 



The last step of our simulation is made calling the PYTHIA software [^g, using as 
input the neutralino decay products above mentioned. 

The signals associated to the first kinematically accessible neutralino, which arise 
from its single R-parity violating production as well as its R-parity conserving pair 
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production are listed in table 1. This table shows the final signals for the x'^t and 
XX production with the subsequent x decay. One of the cleanest and most interest- 
ing signals that can be studied are the events with missing transverse momentum + 
acoplanar muons {"^t + These can be produced through either process a) or 

process b) as shown in table 1. 

The main source of background for this signal is the Z decay to fi~^fi~ with the 
radiation of a 7 which may escape detection. This background has to be eliminated 
through suitable cuts. For definiteness we have imposed the cuts used by the OPAL 
experiment for their search for acoplanar dilepton events [ p^ : 

• we select events with two muons with at least for one of the muons obeying | cos 6'| 
less than 0.7. 

• the energy of each muon has to be greater than a 6% of the beam energy. 

• the missing transverse momentum in the event must exceed 6% of the beam 
energy, jj^t > 3 GeV. 

• the acoplanarity angle (the angle between the projected momenta of the two 
muons in the plane orthogonal to the beam direction) must exceed 20°. 

The OPAL experiment did not find any acoplanar muon pair event passing these 
cuts in the data sample analysed in ref. 



Using our Monte Carlo generator we did the calculation of the kinematic distri- 
butions before and after applying these cuts. This was used to determine the detection 
efficiencies associated to the x^ and XX- They lie in the range ~ 25% — 18% for neu- 
tralino masses in the range ^ 256*6^ — 40GeV for the case of xX production. For 
the single production case {x^) we have found efficiencies in the range 11% — 5% for 
^ 40^61^ - SOGeV. 



7 Analysis and Results 

Given a sample of LEP data collected at the Z peak, the study of the signal described in 
the previous section allows us to determine the corresponding experimental limits on the 
values of the relevant R parity- violating couplings versus neutralino mass. In order to il- 
lustrate the procedure we will use the last data published by the ALEPH collaboration. 
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number hadronic events _i 
corresponding to an integrated luminosity Lint = ; — ; \ — = 82 pb 

(j(e+e^ — > haarons) 

||25|| . This in turn corresponds to 1.94 x 10^ hadronic decays of the Z. Of course this is 



justified only under the assumption that the experimental cuts of the previous section 
are enough to eliminate all relevant background, or that the detection efficiency for our 
signal is not reduced by further cuts that might be needed. Clearly, in order to obtain 
rigorous limits for a given data sample collected in a given experiment, one would have 
to check whether these assumptions are true by means of more detailed simulation 
studies of the corresponding detector features as well as the corresponding background 
for the given luminosity. 



While we await for a more complete statistics to be analysed |^ we find it useful 
to illustrate the sensitivity to the basic parameters of our R-parity violating models 
which has already been achieved with the data samples collected. For such illustrative 
purposes we make use only of the cleanest leptonic signal and use the published ALEPH 
data given in ref. As usual, in order to obtain a 95%CL limit on some parameter, 
we impose the condition 

3 > iVe^pi (38) 

Nexpt being the number of expected events for our signal, when no events of the desired 
type have been observed. 

For the single production process e^e~ — > xz/ with the decay x ~^ z^r/U^/U~ the 
expected number -px + /i+yU~events is given as 

Nexptix^) = c^(e+e" x^)BR{x z^r/U+/i")e^^ Unt (39) 

where e-^y is the detection efficiency, obtained from the R-parity breaking generator 
described before. 

Using the expression for the cross section in eq. ( pT]) we can write 



BR{X l^rfJ'^fJ'')(^xi^ Lint (40) 

where xz = m^/mz- 

In addition, the relation between the coupling O^ls and the BR(Z — * x^) is given 

by 

BRiZ ^ X.) ^ ^O- ^ (l - |4 . i.|) (41) 
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From eq. (^Oj), eq. (^) and eq. ( P5| ) one can obtain a 95%CL limit on the R-parity 
breaking observable BR{Z x^)BR{x J^r/W^/U") as a function of the x mass. This 
is shown in figures 2 and 4. 

Here we should stress that this production mode is characteristic of models with 
spontaneous violation of R-parity [§, ^, or models that parametrize it through an 



effective bilinear superpotential term eiH ||T6|. It is absent in most models of explicitly 



broken R-parity, such as the one considered in the analysis presented by the ALEPH 
collaboration in ref. ||25|| . 



Notice that the R-parity-violating parameter e is directly correlated to the mass 
of the tau neutrino z/^ . As a result it is correspondingly restricted by cosmological 
Big Bang nucleosynthesis ||3^. The corresponding limits are rather stringent ||31| and 



may not allow m,y^ above a few hundred KeV. In this case the signal displayed in Fig. 
4 would never reach 10"^. This cosmological bound can be avoided in models which 
contain neutrino decay or annihilation channels beyond those present within 



the standard model. For this reason the spontaneously broken R-parity models are 
preferred, as they allow the maximum m^^ values permitted by laboratory experiments 
p3| due to the z/^ decays and/or annihilations to majorons. However, as seen in Fig. 



4, the maximum value of our signal is only a bit larger than 10 
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We conclude that in models such as the ones in ref. [T^, ^ one can probe the 
spontaneous violation of R-parity in the single-production mode for ^ 40 GeV. 
However, as already mentioned, these values of the signal rate are hard to reconcile 
with cosmological limits. The corresponding rates for a cosmologically acceptable 
Ur mass lie well bellow the dotted curve in Fig. (^), below 10~^ for all values of the 
neutralino mass, thus too small a rate to be of interest. Including the majoron, whose 
existence is expected in any model where the spontaneous violation of R-parity is 
reahzed in the minimal SU{2) (g) f/(l) gauge structure ||, changes the situation in 
two ways. First, it improves the allowed signal rates because it allows z/^ masses as 
large as present laboratory limits. Unfortunately there is a counter-effect that decreases 
the expected signal rates, because the presence of the invisible channel eq. (Q) tends 
to dilute the x branching ratio into muons. The net result leads to the dotted curve 
in Fig. d). 

The lightest neutralino x ^lay however be light enough to be pair-produced in 
the standard R-parity-conserving process e~^e~ — XX- Even in this case of standard 
production, the violation of R-parity can provide visible signals from the subsequent 



13 



decays of the neutralinos. The first possibihty to consider here is the case where both 
neutrahnos decay visibly, e.g. into, Ur^'^fi'. We do not consider this possibihty as it 
is similar to the one used recently by ALEPH and the corresponding sensitivities may 
be estimated by re-scaling the results of ref. ||25|. Thus we choose concentrate on the 
novel possibility that one of the neutralinos decays to i'rfj^~^fi~ while the other decays 
invisibly, which is more characteristic of models with spontaneous violation of R-parity. 
The number of expected + Z^'*'/^" events in this case is given by 

Nexptixx) = o-(e+e" xx)'^BR{x invisible) z/^/i+/i")e^^ L^nt (42) 

so that from eq. ( P7[ ) we obtain 

2BR{x ^rfJ'^ fJ^—)BR{x ~^ invisible) ej(,^ L^nt (43) 
and the corresponding expression for the Z — > xX branching ratio is 

BRiZ^XX) = i0?„^(l-44)'" (44) 



From these last expressions and eq. (|38D it is possible to obtain an illustrative 



95%CL limit on BR{Z xx)BR{x ~^ ^'r/i^/i )BR{x ~^ invisible) as a function of 
the X mass, as shown in figures 3 and 5. 



8 Discussion 

Using the integrated luminosity corresponding to the last published ALEPH data, L = 
82 pb^^, we have illustrated how the existing data gathered by the LEP collaborations 
at the Z peak are sufficient to start probing in a theoretically meaningful way the mass 
and couplings of the lowest lying neutral supersymmetric fermions in spontaneously 
broken R-parity models. We have determined the corresponding regions of sensitivity, 
illustrated in fig. 2-5, both for the case of single production as well as for the case 
where the lightest neutralino can be pair produced. The theoretical significance of 
these decays is illustrated in Fig. (|^). 

From curve b in this Fig. (^ one sees that the expected signal rate for a neutralino 
of 40 GeV in the model of ref. could be tested if the luminosities of the four 
LEP experiments are added, even if only leptonic + jJ'^ H~ and ]()t + e^e~ channels 
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are considered. This constraint would be significant from tlie point of view of the 
underlying model. Indeed, as can be seen from Fig. (j^), this corresponds to values of 
the mass close to an MeV, for which other processes such as Z — > + |0, 



would be expected to be sizeable. Although the presence of the majoron leads to novel 
decays of standard model particles, e.g. yU — > e+ J and ri+J decays ||TT|, of relevance 



for intense muon source studies at PSI or for a future tau-charm factory as well 
as novel Z decays such as Z — > 7 + J |jT^, its overall effect insofar as the signal rate for 
+ /^"""/^"events arising from production at the Z peak is concerned is to decrease. 
One advantage is that the rate there is no conflict between the signal rate (dotted 
curve in Fig. 3) and the primordial Helium abundance. However, there is clearly a 
wide range of parameters still unconstrained by this single production mechanism. In 
contrast, if kinematically accessible, neutralino pair-production can place more severe 
restrictions, especially on models where the spontaneous violation of R-parity where 
the majoron is absent M or in an effective model where the violation of R-parity is 



introduced explicitly via the bilinear superpotential term iHu The constraints 

derived in this case should be important in relation to searches at higher energies such 
as LEP200, the NLC or the LHC where one expects mostly SUSY pair production to 
yield sizeable event rates. 

Substantial improvements are expected from the use of e~^e~, t^t~ and di-jet plus 
missing momentum event topologies. The results presented here should encourage one 
to perform more detailed and complete background studies and simulations covering 
other R-parity violating signals and improved integrated luminosities already attained 
at LEP From our point of view it would be desirable to have additional runs at 
the Z peak, considering the fact that large areas of parameter space still remain open, 
where the neutrahnos are light enough to be produced in Z decays. 
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A 3-Body Neutralino Decays 



In models with spontaneous R-parity violation, the neutral fermions eigenstates are 
mixed. We denote them collectively by Xi with [i = 1, . . . , 7). The neutrinos corre- 
spond in this notation to i = 1, 2, 3 and the usual neutralinos to i = 4, ... ,7. In the 
charged fermion sector, the leptons are mixed with the charginos of the MSSM. We 
denote them by x~ with {i = 1,...,5). Again, the indices {i = 1,2,3) correspond, 
respectively to e~, /i~ and r^, while the indices i = 4, 5 denote the usual charginos. 
The quarks do not mix with the other fermions. We will consider separately the cases 
where all the particles in the final state have the same charge and the case where they 
can have different charges. 



A.l Same-Charge Final States 



Here the decay we consider is into neutral final states 



(45) 



and proceeds via the neutral current. Due to the Majorana nature of the neutral 
fermions, there are 3 distinct diagrams giving rise to some interference terms. The 
final result for the decay, within the approximation that we neglect all masses of the 
final fermions, is given by 



RJ 



GlMf 
48 7r3 



(Ci + C2 + C3) f{xz) 

+ 2 (Ci2 + Ci3 + C23) g{xz, Xz) 



1 

S p 



(46) 



where 



Cl 
C2 

C3 
C12 

Cl3 
C23 



^ '2 
' Rkm 



'2 
' Rjm 

'2 
'Rim 



(q"2 , q"2 \ (q"2 , q"2 

{pLki + ^Rki) ipLjm + 
(q"2 , q"2 \ (q"2 , q"2 

n" n" n" n" + n" n" n" n" 

^ Lkm^ Lji^ Lki^ Ljm ^ ^ Rkm^ Rji^ Rki^ Rjr 



n" n" n" n" + n" n" n" n" 

^ Lkm^ Rji^ Lkj^ Lim ^ ^ Rkm^ Lji^ Rkj^ Rim 



n" n" n" n" + n" n" n" n" 

^ Lki^ Lkj^ Rjm^ Rim ^ ^ Rki^ Rkj^ Ljm^ Lim 



(47) 



Sp is the symmetry factor for identical particles in the final state, xz = Mi/Mz, and 
the functions f{x) and g{x,y) are given in the Appendix of ref. [^]. The coupling 
matrices are given in section 3. 
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A. 2 Different- Charge Final States 



In this case we consider 3-body final states where two are oppositely charged, that is 

X° ^ + Xk+Xin (48) 

This decay can proceed via both charged and neutral current. Also in the charged 
current case, the Majorana nature of the neutralino implies the existence of two distinct 
diagrams. The final result for the decay, within the approximation that we neglect all 
masses of the final fermions, is given by 



Clf{xz) + (C2 + C3) f{xw) 



GlMf 
48 7r3 



+ 2 (ci2 + C13) g{xz, xw) + 2c23g{xw, xw) 



(49) 



where xw = Mi/Mw, the functions /(x) and g{x,y) are as before, and the coefficients 
Cj are now given by 

ci = (o?. + o2Mo'l^ + o'' 



Lji ' Rji J y^Lkm ' ^ Rkm j 
'^2 = {ylki + ^Rki) {-^Ljm + -^Rjm) 
C3 = iylkj + ^Rkj) i^-^Lim + -^Rim) 
C12 = O'^k^O'^jiYLkiXLjm + O'j^^^O'l^jjYRkiXRjm 
Cl3 = - ip'Lk^O'^j^YLkjXLim + 0'ji^.j^O'[jjYjikjXiiim) 

C23 = yLkiYLkjXnjmXnim + YRkiYRkjX LjmX Hm (50) 



with 



YLij = 


1 

V2 




Yrij = 


1 

V2 
1 

V2 


Kmj 


Xiij = 




Xuij = 


1 

72 


KRji 



(51) 

In all the previous expressions, the coupling matrices are taken to be real. The sign 
factors ej and rji are introduced as explained in section 3. 

Before we close this discussion, we notice that Eq. |^ can also be applied, with obvious 
modifications, to the case of quarks in the final state. For definiteness consider the 
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process 



'X.i 



Xi+u + u 



(52) 



It is clear that in this case we have only the Z exchange diagram. The formula of 
Eq. ^ it is still valid with the following values for the coefficients q 

2 



Cl 
C2 



(03. 



r)"2 



^ - Qu sin^ 9^ 



C3 = C12 = Ci3 = C23 = 



(53) 



Notice that the neutralino decay formulas given above correct those previously given 
in the Appendix of ref. which were not correct for the case of decays involving 
Majorana fermions, such as x — Sur- However the functions f{x) and g{x,y) are the 



same as in ref. p4 |. 
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■T ^ ' V / 


i)T + + (27eis + r) 


T^T + T + 2jets 


h + /+^r + [tI) 


'pT + 2 jets 


i)T + 




i>T + 4jets 




+ 4jeis + T 




jf^T + 2jets + (rZ) 




+ 2 jets 




4jets + rr 




/^T + 2 jets + r 




+ 2 jets + t(t Z) 




;^T(r /)(r 









Table 1: Final signals arising from neutralino pair-production (left column) as well 
single production (right column) . The total charge of the particles between parenthesis 
must be zero. 
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Figure 1: Detection efficiencies for the -pr + l^'^l^ signal associated to ^"^^ XX 
production channels. 
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Figure 2: Region of sensitivity obtained at the 95% C.L. for BR{Z — > x^)BR{x 

as a function of the hghtest neutrahno mass m^o. This is derived from searches 
of ^^T+Zi"^//" events that would arise from single neutralino production at LEP, followed 
by X — > decays. 
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m,o (GeV) 

Figure 3: Region of sensitivity obtained at the 95% C.L. for BR{Z XX){x ~^ 
/i'^ li~i')BR{x — > invisible) as a function of the hghtest neutrahno mass. This is derived 
from searches of + A*"*"/^" events that would arise from neutralino pair production at 
LEP, with one neutralino decaying invisibly and the other decaying as % — > ^'^^"v. 
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m,o (GeV) 

Figure 4: Comparison of the attainable limits on BR[Z — > x^)BR{x ~^ 
versus the lightest ncutralino mass, with the maximum theoretical values expected in 
different R-parity breaking models. The solid line (a) is the same as in Fig. 1, while (b) 
corresponds to the improvement expected from including the e'^e~i' channel, as well 
as the combined statistics of the four LEP experiments. The dashed line corresponds 
to the model of ref. [16] allowing m^^ values as large as the present laboratory limit 
of ref. [23], while the dotted one is calculated in the spontaneous R-parity- violation 
model of ref. [5] . 
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m,o (GeV) 

Figure 5: Comparison of the attainable limits on BR{Z — > xx)BR{x ~^ A*^A*~^) 
versus the lightest ncutralino mass, with the maximum theoretical values expected in 
different R-parity breaking models. The solid line (a) is the same as in Fig. 2, while (b) 
corresponds to the improvement expected from including the e+e"!/ channel, as well 
as the combined statistics of the four LEP experiments. The dashed line corresponds 
to the model in ref. [16] allowing m^^ values as large as the present limit, the dotted 
one does implement the restriction on rrii,^ suggested by nucleosynthesis, and the dash- 
dotted one is calculated in the model of ref. [5] . 
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m,,(MeV) 

Figure 6: Attainable values for the R-parity breaking strength BR{Z x^) versus 
the T-neutrino mass in the model of ref. [5]. 
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